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ABSTRACT 


Dynamically active planetary systems orbit a significant fraction of white dwarf stars. These stars 
often exhibit surface metals accreted from debris disks, which are detected through infrared excess 
or transiting structures. However, the full journey of a planetesimal from star-grazing orbit to final 
dissolution in the host star is poorly understood. Here, we report the discovery that the cool metal 
polluted star WD 0816-310 has cannibalized heavy elements from a planetary body similar in size to 
Vesta, and where accretion and horizontal mixing processes have clearly been controlled by the stellar 
magnetic field. Our observations unveil periodic and synchronized variations in metal line strength 
and magnetic field intensity, implying a correlation between the local surface density of metals and the 
magnetic field structure. Specifically, the data point to a likely persistent concentration of metals near 
a magnetic pole. These findings demonstrate that magnetic fields may play a fundamental role in the 
final stages of exoplanetary bodies that are recycled into their white dwarf hosts. 
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1. INTRODUCTION 


White dwarf stars exhibit multiple hallmarks of their 
remnant yet vibrant planetary systems, including cir- 
cumstellar debris that sometimes transits the host (Far- 
ihi 2016; Guidry et al. 2021). As the planetary debris 
is accreted onto the white dwarf, it results in a tempo- 
rary phase during which photospheric metals can be de- 
tected, often referred to as pollution (Zuckerman et al. 
2003). Because of the downward settling of heavy el- 
ements in their high-gravity atmospheres, white dwarf 
spectra normally exhibit only hydrogen or helium, but 
a significant fraction also show traces of these accreted 
metals (Zuckerman et al. 2010; Koester et al. 2014). 
This phenomenon is now well understood to be the result 
of accretion of planetary material, where the relative el- 
emental abundances measured in polluted white dwarfs 
provide a gateway to studies of extrasolar geochemistry 
and planet formation (Klein et al. 2021; Bonsor et al. 
2023). 

The observed debris disks are attributed to star- 
grazing planetesimals that are tidally fragmented, but 
their evolution is understood only in the broadest terms 
(Malamud & Perets 2020; Brouwers et al. 2022). For 


example, the necessary drastic compaction of the semi- 
major axis, from several astronomical units to periastra 
within a solar radius, is an active area of research in 
which collisions are likely to play a central role (Veras 
et al. 2015; Malamud et al. 2021). In some cases mag- 
netic fields may assist in orbital circularization for some 
disk constituents (Hogg et al. 2021; Zhang et al. 2021). 
The accretion of disk particles onto the star is likely to 
be driven by Poynting-Robertson drag, and possibly en- 
hanced by solid-gas coupling or the collisional destruc- 
tion of solids (Rafikov 2011; Kenyon & Bromley 2017). 
In this process, magnetic fields may play a fundamental 
role by re-directing the gaseous accretion stream (Met- 
zger et al. 2012; Farihi et al. 2018). 

In the local 20pc volume of white dwarfs, at least 
15 per cent are metal-polluted (Hollands et al. 2018), 
while over 20 per cent host a detectable magnetic field 
(Bagnulo & Landstreet 2021), with a significant overlap 
between these two populations (Farihi et al. 2011; Hol- 
lands et al. 2018; Bagnulo & Landstreet 2019a; Kawka 
et al. 2019). It is well known that in main-sequence 
stars, the presence of a globally organised magnetic field 
enables and preserves a non-homogeneous distribution 
of the chemical elements of the photosphere (e.g. Do- 
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nati & Landstreet 2009). An analogous, but physically 
distinct phenomenon has been observed in three mag- 
netic white dwarfs that have both hydrogen and helium 
lines in their spectra (Achilleos et al. 1992; Wesemael 
et al. 2001; Moss et al. 2024), but, perhaps surprisingly, 
it has hardly been investigated in metal-polluted stars. 

WD 0816-310 is a helium atmosphere white dwarf at 
19.4pc distance, exhibiting a DZ spectral type with 
strong absorption lines from multiple metal species, 
most prominently Ca, Mg, and Fe (Subasavage et al. 
2008). Spectropolarimetric observations have detected 
a variable magnetic field (Bagnulo & Landstreet 2019a). 
A spectroscopic study by Kawka et al. (2021) argued 
that the metal abundances in WD 0816-310 had changed 
between two observations taken 10 yr apart, and that 
this change was probably a consequence of the surface 
magnetic field. Here, we report the results of a spec- 
tropolarimetric monitoring of this star, aimed at discov- 
ering whether there is a relationship between the struc- 
ture of the magnetic field and the distribution of the 
metal pollutants on the surface. 


2. OBSERVATIONS 


New circular spectropolarimetric observations of 
WD 0816-310 were obtained in 2023 February and 
March with the FORS2 instrument (Appenzeller et al. 
1998) of the European Southern Observatory Very 
Large Telescope at Cerro Paranal. FORS2 is a 
multi-purpose instrument capable of imaging and spec- 
troscopy, equipped with polarimetric optics. Five obser- 
vations were obtained on four different nights with the 
grism 1200B, covering the spectral range 3700-5100 A, 
at resolving power R ~ 1400. Data were used to deter- 
mine the longitudinal component of the magnetic field 
averaged over the stellar disk, or mean longitudinal field 
(Bz), and to search for line strength variability (see Ap- 
pendix A). The dataset shows large changes in the value 
of both (B,) and the apparent strength of a number of 
absorption lines over an interval of a few days (see Fig- 
ure 1). In particular, the variable lines are wider and 
deeper when (B+) is large and positive (2023 February 
24 and 25), and weaker when (B-,) is close to zero (2023 
February 27 and March 26). The magnetic field mea- 
surements, observing dates, and instrument setup are 
provided in Table 1. 


3. RESULTS 
3.1. Rotation period and magnetic field morphology 


The measured longitudinal field variability is naturally 
explained in terms of a static field that is not symmetric 
about the stellar rotation axis, so that the observer sees 
a varying field as the star rotates. The simplest possible 


magnetic configuration is a dipole field inclined at an 
angle 6 to the white dwarf rotation axis, which in turn 
is inclined at an angle i to the line of sight. Accord- 
ing to this oblique rotator model, originally introduced 
to explain the field variations observed in the chemically 
peculiar stars of the upper main sequence (Stibbs 1950), 
the longitudinal field plotted as a function of rotational 
phase has a simple sinusoidal shape. In fact, contribu- 
tions to (B+) from higher order multipoles tend to can- 
cel out when averaged over the stellar disk, so that (Bz) 
measurements are mainly sensitive to the contribution 
from the dipolar component (Bagnulo et al. 1996), if it 
exists. 

Sinusoidal curves with several periods between 2 and 
16d fit the (B+) values satisfactorily. While the data 
cannot uniquely determine the spin period, they never- 
theless demonstrate that during the rotational cycle, the 
positive magnetic pole and the magnetic equator cross 
the line of sight, and at certain rotational phases, the 
negative magnetic pole transits the visible stellar hemi- 
sphere. In spite of a periodogram with many candidate 
signals (and aliases of each), a good fit to the vary- 
ing longitudinal field measurements, assuming a dipolar 
field, depends only weakly on the choice of the rotational 
period. Panel (a) of Figure 2 shows the best fit to the 
(B.) values obtained by adopting a rotation period of 
10.89426d. A plausible and simple model consists of a 
star with i ~ 6 ~ 70°, and dipolar field strength at the 
pole ~ 140kG (see Appendix B). Panel (b) shows the 
equivalent width (EW in Table 1) of the Fe1 lines in the 
region 4357 — 4425 A plotted with the same period and 
ephemeris. Only one spectrum probes the region charac- 
terized by a negative field. Additional data are required 
to find the exact rotation period, to improve the mod- 
elling, and to establish how the equivalent widths and 
(Bz) correlate in the region characterized by a negative 
longitudinal field. 

There are five available sectors of PDCSAP photome- 
try available for WD 0816-310 from the Transiting Ex- 
oplanet Survey Satellite, (TESS) spanning close to 4 yr 
between Sectors 07 and 61. However, these data reveal 
no significant periodogram signals above an amplitude 
of 0.58 per cent, and none that convincingly stand out 
from adjacent signal peaks. Furthermore, only 10-15 per 
cent of the PDCSAP flux is from the white dwarf, while 
over 85 per cent is from other sources. Thus, no rotation 
period can be confidently extracted from these data. 


3.2. Metal line strength variability 


A comparison between a FORS1 spectrum obtained in 
2007 and an X-Shooter spectrum obtained in 2017 re- 
vealed variations in the strength of the lines correspond- 
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Figure 1. Stokes spectra of WD 0816-310, with flux normalized to the continuum shown in panel (a) and V/I plotted in 
panel (b). The unsaturated lines are stronger when the longitudinal field is near maximum, and weaker when the field is close 
to zero. The presence of a field is revealed by the non-zero Stokes V profile, as detailed in Appendix A. For visual clarity, in 
panel (b) only the V/I profiles obtained on 2023 February 24 and March 26 are shown, and their uncertainties are plotted in 
cyan, and offset by —7.5%. Line identification from Kawka et al. (2021). 


Table 1. Magnetic field measurements of WD 0816—310 obtained from the (saturated) Ca 1 lines and EW from the lines in 
the range 4357-4425 A. 


INSTR. Grism/ DATE UT Exp (Bz) (kG) EW (A) Ref. 
Grating yyyy-mm-dd hh:mm (s) (Caml) (Fe 1) 

FORS1 600B 2007-10-22 08:40 960 —27.0+3.0 3.35+0.1 BLI19 

ISIS R600B 2019-04-19 21:16 3600 43.8+15.0 5.5440.4 BLI9 


FORS2 1200B 2023-02-24 05:48 2700 47.9+3.9 5.60+0.1 this work 
FORS2 1200B 2023-02-25 02:20 2700 42.2+3.0  5.14+0.1 this work 
FORS2 1200B 2023-02-25 05:12 2700 42.642.6  5.04+0.1 this work 
FORS2 1200B 2023-02-27 05:14 2700 —1.4+2.6 3.27+0.1 this work 
FORS2 1200B 2023-03-26 00:34 2700 1.642.4 3.17+0.1 this work 


References—BL19: Bagnulo & Landstreet (2019a). 
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Figure 2. Panel (a): the best fit sinusoidal model to the 
(Bz) data of Table 1, having adopted a rotation period of 
10.89426 d. The zero point of the rotational phase is set ar- 
bitrarily at JD=2460000.0. Other values of the rotational 
period in the range between 2-16d lead to equally satisfac- 
tory fits with longitudinal curves of similar amplitude. Panel 
(b): the equivalent width of the lines in the region between 
4357 and 4425A phased to the same rotation period and 
ephemeris as (Bz). 
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Figure 3. All VLT flux spectra of WD 0816-310 analyzed 
in this paper, normalized to the continuum, in a wavelength 
range dominated by Cat and Fe" lines. Note that two spec- 
tra obtained on 2023 February 25, nearly overlapping, are 
plotted using the same color. 


ing to Na, Mg, Ca, Cr, Mn, Fe, and Ni (Kawka et al. 
2021). Atmospheric modeling demonstrated the elemen- 
tal abundance changes were clustered around a 0.3 dex 
offset between the lower values estimated from the 2007 
FORS1 spectrum, and the higher values found with the 
X-Shooter data taken 10 yr later (Kawka et al. 2021). 
The essential conclusion drawn from our new data is that 
these line strength changes occur on a timescale of days, 
are strongly correlated with the longitudinal component 
of the magnetic field, and therefore with the rotation of 
the star. By ruling out the possibility that the abun- 
dance changes might be caused by recent accretion (see 
also Sect. 3.3), these data represent the first clear de- 
tection of abundance spots on a cool polluted magnetic 
white dwarf. 

The FORS2 spectra obtained in 2023 February 24 and 
25 show spectral lines of similar strength as those ob- 
served with X-Shooter in 2017; the other two recent 
FORS2 spectra, for which we measure a (B,) close 
to zero (2023 February 27 and 2023 March 26), have 
weaker lines that nearly coincide with those of the ear- 
liest FORS1 spectrum (2007), from which we measured 
a negative (B,) value. It is clear that the two spec- 
tra previously analyzed by Kawka et al. (2021) show 
roughly the extrema of line strengths for the full avail- 
able dataset, and therefore define approximately the 
range of (hemispherically averaged) abundances at the 
surface of WD 0816-310 as seen by the observer as the 
star rotates. All flux spectra in the region 4150-4490 A 
are displayed in Figure 3. 

Local atmospheric abundance variations (or abun- 
dance patches) are commonly observed in magnetic Ap 
stars, and one may be tempted to assume that the same 
mechanism operates in metal polluted white dwarfs, but 
the physics is entirely distinct. At the surface of an Ap 
star, abundance patches result from selective radiative 
levitation of specific ions, guided by the magnetic field, 
causing large local over-abundances of high-Z elements. 
In sufficiently cool white dwarfs, radiative forces are no 
match for gravity, and all heavy elements sink from the 
atmosphere in a short time. Metals observed in the at- 
mospheres of white dwarfs are understood to be the re- 
sult of a relatively recent accretion of material derived 
from asteroid-like bodies injected onto star-grazing or- 
bits, following the gravitational perturbation by a planet 
(e.g. Bonsor et al. 2011). The asteroid is then disrupted 
by tidal forces when it passes within ~ 1 Rọ from the 
star, and the fragments gradually form a circumstellar 
debris disk whose material is accreted (Jura 2003). In 
Tes X 11000K, non-magnetic white dwarfs, 3D atmo- 
spheric models predict that metals will have a homo- 
geneous distribution over the stellar surface (Cunning- 


ham et al. 2021). Before discussing how the observed 
abundance inhomogeneities may have been formed and 
preserved, we examine the accretion history of the star 
as traced by the multiple heavy elements in the atmo- 
sphere. 


3.3. Evidence for a post-accretion epoch 


The photospheric metal-to-metal ratios reveal that 
Mg is notably over-abundant, relative to all heavier ele- 
ments (Ca, Cr, Mn, Fe, Ni), by a factor of 2-4 as com- 
pared to the ratios found in chondritic meteorites, fol- 
lowing a trend that correlates with atomic mass and 
hence diffusion time (Appendix C; Lodders 2003). This 
is a well-understood pattern arising from the chemical 
separation of heavy element species according to their 
distinct sinking timescales. In the decreasing abundance 
phase of an accretion episode, infalling material is ex- 
hausted, and heavier elements such as Fe sink more 
rapidly, while the lighter element Mg remains in the 
photosphere longer (Swan et al. 2019). Even under the 
assumption of steady-state accretion, the Mg overabun- 
dance persists, but there is no geochemical pathway to 
enhance this element over all heavier metals. 

In some polluted white dwarfs, the intrinsic abun- 
dances of the debris exhibit departures from chondrites, 
and in particular can show evidence for differentiated 
material such as observed in planetary crusts and cores 
(Zuckerman et al. 2011; Gansicke et al. 2012; Jura et al. 
2013). However, to the limits of the uncertainties, most 
white dwarfs have accreted planetary bodies that are 
broadly consistent with chondritic, undifferentiated ma- 
terial (Jura & Young 2014; Xu et al. 2019; Doyle et al. 


2023). 
To model the effects of sinking on the observed heavy 
elements, we adopt stellar parameters Tig = 6250K 


and logg = 8.25 (cgs), together with the linear aver- 
age of two sets of metal abundances for Mg, Ca, Cr, 
Mn, Fe, and Ni calculated by Kawka et al. (2021). 
These parameters were then used to define the corre- 
sponding atmospheric diffusion models for helium at- 
mosphere white dwarfs with metals (including overshoot 
and [H/He] = —5; Koester et al. 2020). 

Figure 4 plots the three mass ratios Ca/Mg, Fe/Mg, 
and Fe/Ca corresponding to the parent body of the ac- 
creted material, over time as the three elements sink ac- 
cording to their individual rates. The plot demonstrates 
clearly that all three ratios are generally consistent with 
chondritic values, if accretion ceased between 0.4 and 1.2 
Mg sinking timescales (TMg) in the past, depending on 
the prior accretion history. The best-fitting model cor- 
responds to 0.31 Myr post-accretion, following a steady 
state. 
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Figure 4. Time-reversed, post-accretion evolution of the 
mass ratios Fe/Ca, Ca/Mg, and Fe/Mg for the accreted 
source material in WD 0816-310, normalized to the values of 
chondritic meteorites (Lodders 2003). Time on the x—axis 
corresponds to t—to in Equation C9, increases toward earlier 
epochs, and is normalized by the Mg sinking time. Two cases 
are plotted, where the photospheric mass ratios at to are the 
result of 1) an increasing phase, shown by dashed lines, or 
2) steady-state accretion, plotted as solid lines. All three 
ratios exhibit an under-abundance of the heavier element at 
the current epoch, but as time moves into the past, these ap- 
proach and eventually broadly match those of the primitive 
inner solar system. Comparing the two plotted models, the 
solid lines span a narrower range of times when crossing over 
with the chondritic model, thus supporting a steady state 
of accretion prior to its cessation. An additional figure and 
details are provided in Appendix C. 


Importantly, these three elements behave distinctly in 
the case of planetary differentiation (e.g. by melting): 
Ca is strongly lithophilic and greatly enhanced in plan- 
etary crusts, Fe is strongly siderophilic and dominates 
cores, while Mg is only mildly lithophilic. These dis- 
parate geochemical affinities would readily be revealed 
in accreted material derived from a differentiated par- 
ent body. As none of these tendencies are observed, the 
abundance modeling indicates that WD 0816-310 is pol- 
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luted by material similar to chondritic meteorites, and 
that accretion has ceased. 

The above analysis is self-consistent and agrees with 
prior work on polluted white dwarfs in the decreasing 
abundance phase of accretion (e.g. Swan et al. 2023). 
On their own, these calculations argue that the mag- 
netic field of WD 0816-310 has no significant impact on 
the relative sinking timescales for multiple heavy ele- 
ments. This should not be too surprising, as the ratios 
of sinking timescales for pairs of metals change little over 
wide ranges of effective temperature and surface grav- 
ity (Koester 2009), and thus the results of this anal- 
ysis are robust in the case that magnetism alters the 
absolute sinking times. However, even this is unlikely, 
because atmospheric modeling predicts the field is too 
weak to significantly impact the size of the convection 
zone (Bédard et al. 2017; Cunningham et al. 2021). 

The most robust conclusion possible is that the pho- 
tospheric metals have had at least 0.4 ne or 0.3 Myr 
to diffuse downward as well as horizontally. Further- 
more, the total mass of metals in the convection zone 
requires that accretion had been ongoing for a signifi- 
cant period prior to ending. A minimum mass for the 
disrupted and accreted parent body can be calculated 
by assuming the photospheric Fe is 20 per cent of the 
total mass as in chondrites, resulting in 2 x 1071 g re- 
siding in the star at present. If this material has been 
constantly sinking for roughly one diffusion time, then 
the minimum mass becomes D x 107! g. However, if the 
star experienced steady-state accretion prior to ceasing, 
as suggested by the analysis here, then the mass of the 
parent body can easily be one to two orders of magni- 
tude larger, approaching or exceeding the mass of the 
asteroid Vesta. Future improvements in diffusion mod- 
els will hopefully provide better constraints on accretion 
history (e.g. Heinonen et al. 2020). 


4. THE ORIGIN OF THE METAL PATCHES 


Spectral variations owing to changes of the local abun- 
dance of an element can occur only if the atmosphere 
is composed of at least two elements. Among white 
dwarfs, it is quite rare that both hydrogen and helium 
are comparably abundant and simultaneously exhibit 
absorption features (e.g. spectral type DBA), but in the 
case of external pollution, there are metals and either 
hydrogen or helium (or both). Prior to the study of 
WD 0816-310, variability of the H/He ratio had been 
found only in a few DBA white dwarfs. Among these, 
Feige 7, PG0853+164 and SDSS J091016.48+210554.2 
are the only stars that also show evidence for a (variable) 
magnetic field (Achilleos et al. 1992; Wesemael et al. 
2001; Moss et al. 2024). A few other DBA stars, such as 


G104-27 (Kidder et al. 1992), PG 1210+533 (Bergeron 
et al. 1994), WD 0209+085 (Heber et al. 1997), GD 323 
(Pereira et al. 2005), and ZTF J203349.8+322901.1 (Ca- 
iazzo et al. 2023), show periodic, or at least recurrent, 
variations in their H/He abundance ratios. In these 
stars, the presence of a weak magnetic field has often 
been suggested, but no sensitive magnetic measurements 
have been made. In addition to these cases, the presence 
of undetected helium in the atmosphere of the strongly 
magnetic, cool DAH star G183-35 has been proposed 
(Kilic et al. 2019) to account for the observed variability 
of the hydrogen line strength, caused by a non-uniform 
distribution of hydrogen. 

Feige 7 and SDSS J091016.43+210554.2 are the only 
white dwarfs for which modelling definitively confirms 
the direct influence of a magnetic field on localized vari- 
ations in the atmospheric H/He ratio (Achilleos et al. 
1992; Moss et al. 2024). It is thought that in these hot 
white dwarfs, the magnetic field completely suppresses 
the expected surface convection, and thus prevents hor- 
izontal mixing, thereby enabling the persistence of dis- 
tinct patches with different H/He ratios. However, the 
mechanism that has led to the formation of the observed 
patches is not understood. 

WD 0816-310 is the first well-established analogue to 
Feige 7, and is the first object found to show abundance 
patches among cool, metal-rich white dwarfs. In this 
DZ star, the non-uniformly distributed metals are the 
remnants of a disrupted asteroid or analogous planetary 
body, and were deposited from a circumstellar reservoir 
onto the surface of the white dwarf. We now discuss how 
metal abundance patches may arise from this specific 
situation. 

Sufficiently close to the white dwarf, any circumstellar 
planetary debris will be entirely gaseous, as the equilib- 
rium conditions will eventually exceed the condensation 
temperatures for highly refractory minerals (Rafikov & 
Garmilla 2012; Okuya et al. 2023). In the absence of a 
magnetic field, the accreting gas would be expected to 
arrive on the stellar surface in a band defined by the rela- 
tive orientation of the disk and the white dwarf rotation 
axis. The material would then be expected to be spread 
more or less uniformly over the full white dwarf surface 
by atmospheric and envelope turbulent convection, on 
a timescale much shorter than for vertical diffusion (see 
fig. 11 in Cunningham et al. 2021). 

The presence of the stellar magnetic field may alter 
this picture. According to 3D atmospheric models of 
vertical and horizontal mixing in white dwarfs (Cun- 
ningham et al. 2021, Figure 13), a 100kG field is insuf- 
ficient to inhibit surface convection in a Teg zs 6250K, 
helium atmosphere star such as WD 0816-310. However, 


in contrast to these models, the observations show that 
the ~% 100kG field appears sufficiently strong to inhibit 
horizontal spreading and preserve the non-uniform ar- 
rival pattern of metals, for at least 0.3 Myr (based on 
Figures 4 and 6, Section 3.3 and the Appendix D). 

The mechanism responsible for the initial formation 
of the abundance patches can be tied to the accretion 
process as follows. It is well established that even a 
weak field will truncate a sufficiently-ionized gas disk 
and guide the accretion stream to preferentially reach 
the surface in two patches, one at each magnetic pole 
(Metzger et al. 2012; Farihi et al. 2018). However, based 
on the Teg + 6250 K of WD 0816-310, it is not obvious 
that the inward drifting gas will be sufficiently ionized 
for magnetic control of the accretion flow. Almost all 
of the accreting metal atoms have ionization thresholds 
below 2800 A, where models of cool DZ stars predict flux 
less than 3 per cent of that near 4000 A. Therefore, it is 
unclear if sufficient ultraviolet flux is available to signif- 
icantly ionize the accretion flow. However, we propose 
that a powerful ionization multiplier mechanism will be 
at work in the inner gaseous accretion disk of WD 0816- 
310. 

Because of near conservation of angular momentum, 
in the region between the star and disk where solids 
rapidly sublimate, the circumstellar neutral gas would 
orbit the white dwarf with a velocity between around 
1500 and 3500kms~'. In other words, metal atoms in 
orbit carry typical kinetic energies of the order of 1 MeV. 
As soon as one of these metal atoms becomes ionized, 
it is immediately affected by the magnetic field of the 
star, which ranges in strength from roughly 1kG where 
the disk transitions from solids to gas, to 100kG near 
the star. As a result of the Lorentz force, the ion will 
immediately enter a local circular or spiral orbit, with 
a radius of order meters, around a local field line. It 
will no longer move with the neutral swarm of atoms, 
but will at first continue to move at co-rotation speed 
but in rapidly changing directions. Consequently it will 
collide with neighbouring neutral streaming particles at 
high velocity and energy. These collisions will hardly 
ever restore the rogue ion to neutrality; instead the ion 
will probably achieve higher ionization states. Further- 
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more, most of the collisions will result in ionization of the 
target neutral atoms. The energy available for further 
ionizations from creation of a single seed ion is enough 
to ionize thousands of neutral atoms. Many of these 
new ions will in turn produce still further ions (see Ap- 
pendix D for additional detail). 

In summary, we propose that, even if neither the avail- 
able ultraviolet flux, nor the gas temperature are suffi- 
cient to adequately ionize the disk particles, the pres- 
ence of a tiny number of atoms ionized by the stellar 
flux will progressively lead to a substantial level of ion- 
ization within the gas disk. Consequently, the accreting 
stream will be forced to follow the field lines towards 
the magnetic poles on the surface of the white dwarf. 
In the polar regions, the inhibition of horizontal mixing 
will then facilitate the accumulation of metals. 

It may be surprising that abundance patches in mag- 
netic, metal-polluted white dwarfs have not been discov- 
ered until recently. One reason for this may be the gen- 
eral tendency to obtain only a single spectrum of each 
target white dwarf, or a lack of detailed comparison be- 
tween multiple spectra, where available. Furthermore, 
the presence of abundance patches may well escape de- 
tection unless spectra are obtained with high signal-to- 
noise, e.g. S/N > 100 in the continuum, a level routinely 
achieved when measuring polarization in spectral lines, 
but rarely in large spectroscopic surveys. We predict 
that by adopting a tailored observing strategy, line vari- 
ability will become a frequently observed phenomenon 
in magnetic polluted white dwarfs. 
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APPENDIX 


A. OBSERVATIONS, DATA, AND MAGNETIC FIELD MEASUREMENTS 


The polarization spectra were obtained using the beam swapping technique (Bagnulo et al. 2009), and reduced using 
the FORS2 pipeline (Izzo et al. 2010), although the extraction of the wavelength calibrated beams was made with 
IRAF, and beam recombination was performed using dedicated FORTRAN routines. Intensity (Stokes I) spectra were 
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then divided by the spectrum of a featureless DC white dwarf obtained with the identical instrumental setup, and 
then normalized to a convenient continuum. Owing to the beam-swapping technique, the measurements of the fraction 
of circular polarization (Stokes V divided by the intensity I) are self-calibrated (Bagnulo et al. 2009). The sign of 
Stokes V is defined as in previous works (Bagnulo & Landstreet 2019b). 

The observed fields are sufficiently weak that the profiles of absorption lines in the (R ~ 1000) spectra are virtually 
unchanged from those expected in the absence of a magnetic field. Thus, the mean longitudinal field (B,) is determined 
via a least-squares technique using the correlation between polarization and line slope (Bagnulo et al. 2002) 


21 dI 


= — Jeff Cz À TD (Bz) (A1) 


I 
with 


fs ~ 4.67 x 1078 A“! Go} (A2) 


SC" Aarm, C7 
where e is the electron charge, Mme the electron mass, c the speed of light, and geg the effective Landé factor. We 
note that the derived mean longitudinal field is weighted over the observed stellar hemisphere by the local brightness 
(which decreases towards the stellar limb), and by the local strength of the spectral line(s) used for the measurement. 

In this work, we analyse the newly obtained FORS2 data, and re-visit observations obtained with FORS2 and the 
ISIS instrument of the William Herschel Telescope, which were published in a previous discovery paper (Bagnulo 
& Landstreet 2019a). All of the observing series of sub-exposures were obtained by setting the retarder waveplate 
at sequential position angles —45,+45,+45,—45°. Because white dwarfs may have short rotation periods, we also 
split each observing series into two pairs of sub-exposures, to search for short-term variability. Because no evidence 
of variation within a 1h observing series was found, the Stokes profiles obtained from each complete series of four 
sub-exposures were used in this work. 

The method of determining (B+) through Eq. Al has been well tested on many kinds of stars, including white dwarfs 
(Bagnulo et al. 2012; Bagnulo & Landstreet 2018). Overall, it leads to (B,) values broadly in agreement with those 
obtained with other methods, but it is invariably found that results depend on the specific choice of spectral lines or 
instrument setup (Landstreet et al. 2014). In the case of WD 0816-310, it was found that (B+) values obtained from 
the highly saturated Cat lines and (B,) values obtained from the weaker, variable lines around 4214-4415 A show 
some systematic differences around the epoch of maximum field. For example, for the data obtained on 2023 February 
24, the analysis of Ca 1 lines yields (B,) = +48+4kG, while the analysis of the spectral lines in the range 4214-4415 A 
yields (Bz) = +58 + 3kG. 

These differences may arise in part from simplifications inherent in the correlation method, which is correct when 
applied to single lines, but is only approximately correct for a spectral region crowded with numerous blended lines. 
Furthermore, the Call lines are strongly saturated, and their shape changes little as the ion abundance varies. While 
is it not known how the strong saturation of the two Ca II lines weight the local contributions to the (B,) measurement 
over the visible hemisphere, the similarity of the local Cam line profiles in the observations obtained at different 
rotation phases may support a roughly uniform weight. The weaker lines around 4214-4415 A are instead sensitive 
to small variations of the relevant abundance. A higher weight will be given to the values of the magnetic fields that 
characterizes the regions where the elements are more abundant. As metal elements appear more abundant around 
the north magnetic poles, they probe with higher weight the regions of the stellar hemisphere characterized by the 
highest field strength. 

Of course, a significantly better characterization of the magnetic field may be achieved with proper modelling of 
the observed Stokes parameters, in a way similar to what is currently done for Ap/Bp stars. Such modelling is not 
possible with the limited set of low-resolution polarization spectra presently available. 

In Table 1, we report the field estimates obtained from the combined spectral region 3906-4007 A that contains the 
strong Call lines, using ge = 1.25. Figure 5 shows the correlation of Eq. Al for two observing datasets. These are the 
values adopted in the paper, keeping in mind that they are field measurements based on a particular spatial sampling 
as described above. 


B. AN APPROXIMATE MODEL OF THE STELLAR MAGNETIC FIELD 


The simplest possible model of the stellar magnetic field is a dipole of polar field strength Bp inclined at an obliquity 
angle 6 to the white dwarf rotation axis, which in turn is inclined at an angle 7 to the line of sight. Both angles range 
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Figure 5. Eq. Al correlations for the observing dates and spectral ranges as given in the figure panel labels. The slopes of 
the fitted lines are proportional to the mean longitudinal magnetic field. 


from 0 to 180° (Landolfi et al. 1998), and the variation of the longitudinal field with phase has a simple sinusoidal 
shape (Stibbs 1950) 
(Bz) ~ 0.31B, | cosicos 8 + sin isin B cos(f — bail (B3) 


where f is the rotation phase and fo a zero-phase point such that (B,) is maximum when f = fo. Equation B3 shows 
that if the dipolar field is a reasonable approximation, then the magnetic curve may be fit by a first order Fourier 
expansion, an assumption that is adopted throughout this work. 

The best fit to the (B,) data of Table 1 is then obtained with a rotation period of 10.89426d. This value is not 
unique, and would be different if (B,) were obtained from a different interval range. Nevertheless, given the fact that 
the longitudinal field values has been close to its maximum for about 2 days, and dropped to zero after another 2 days, 
it is likely that the rotation period is between a few days and ~ 2 weeks. 

The (B,) curve of Eq. B3 is degenerate, and does not change either if i and 6 values are swapped, nor if i and 8 
transform into 180° —i and 180° — 0. respectively. The latter transformations correspond to a reversal of the rotation 
direction; therefore, for a physical understanding of the model, it is sufficient to consider only 0° < i < 180° and 
0° < 8 < 90°. A feature of the limited dataset is that (B,)™** > 0 > (B.)™™ and |(B.)™**| > |(B.)™™|. This implies 
that i + 8 must exceed 90° in order for the line of sight to intersect both magnetic field hemispheres, and that 7 < 90° 
(see fig. 1c of Landolfi et al. 1997). We also note that our data allow a rough estimate of the (B,) curve extrema, 
which are reached, respectively, at f — fo = 0° ((Bz) ~ 48kG), and at f — fo = 180°0 ((B,) ~ —39 kG). We define 

„= (Bam _ conti + 8) T 

(B,)™*  cos(i— 8) 
(using Eq. Al), and obtain r ~ —0.8. There are an infinite number of pairs (7,3) satisfying this condition. However, 
some constraints come from the lack of Zeeman splitting in FORS2 data, and from the estimate of the mean field 
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Figure 6. The photospheric mass ratios Z/Mg for WD 0816-310, plotted as filled red circles with error bars, ordered by atomic 
mass, and compared to models. The ratios are the linear average of the two sets of elemental abundances determined from 
atmospheric model fitting (table 3; Kawka et al. 2021), where conservative errors of 0.1dex have been adopted based on the 
internal consistency of the published values. The data are plotted relative to CI chondrites, which are frequently identified as a 
good match for the parent bodies of white dwarf pollution (see main text). There are four models shown for comparison, with 
varying degrees of agreement, all based on the accretion of chondritic material. Model 1 ()open black stars) is simply chondritic 
Z/Mg, such as might be directly observed in the star during an increasing abundance phase. Model 2 (open pink circles) is 
similar, but where these ratios are modified by steady state accretion. Model 3 (open blue squares) assumes a decreasing phase 
of 0.76 Myr follows an increasing phase, and model 4 (filled blue squares) represents a decreasing phase of 0.31 Myr after a 
steady state. 


modulus (|B|) obtained from X-Shooter (Kawka et al. 2021), where (DI ~ 92kG. A reasonable approximation of the 
(|B|) curve for a dipole field is (Hensberge et al. 1977) 


(Bl) = 0.64B, (1 + 0.25 [cosicos 8 + sinisin B cos( f — fo)|’) (B5) 


Therefore, the X-Shooter measured (|B|) value provides the constraint 115kG < Bp < 145kG. Solutions that keep 
alos the maximum of the longitundinal field > 40kG are those with i = 70° + 10° and 6 = 70° + 10°, for example 
i ~ p ~ 70° and a polar field strength Bp ~ 130 — 145 kG. 

A careful comparison of the core shapes of the Catt H and K lines in the FORS 1200B spectra with the X-Shooter 
spectra indicates that none of our spectra are consistent with (| B|} > 150kG. This result suggests that by expanding 
the field structure in a multipole series, the dipole component would still dominate, as the model dipole component 
accounts for about 100 kG of the observed H and K line core broadening. 
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C. ACCRETION-DIFFUSION CALCULATIONS FOR HEAVY ELEMENTS 


The calculations performed in this paper follow the formal treatment of debris disk accretion and diffusion theory 
(Koester 2009; Jura et al. 2009). Specifically, there are three distinct phases that are considered for time t and 
representative sinking timescale 7; an increasing abundance phase when accretion begins and t < 7, followed by a 
steady-state where accretion is ongoing, and abundances are constant for t > 7, and finally a decreasing abundance 
phase after accretion has ended (see e.g. fig. 2 of Koester 2009). 

There are accurate approximations to the full set of differential equations, and below it is assumed that accretion 
rate is a constant, and the transition to the decreasing phase is abrupt. The instantaneous accretion rate for element 
z in the steady state is given by 


. X, M, 
M, = e, (C6) 
Tz 
where 
Mz, 
X, = — x 10/2/"el (C7) 
MHe 


is the mass fraction of the element z within Mevz, the mass of the convection zone or fully mixed layer of the star, and 
[Z/He] is its logarithmic number abundance relative to helium. As detailed in the main paper, the stellar properties 
(Mevz,Tz) were taken from published model grids, and [Z/He] from published atmospheric model fitting of WD 0816- 
310 (Kawka et al. 2021). 

In the increasing abundance phase, the mass ratio for any pair of elements in the star, X,1/Xz2 (e.g. Ca/Fe) 
is identical to that within the parent material, as diffusion has yet to play a role. Once a steady-state has been 
established, this ratio is now modified as follows: 


Mz = Xz1 HES (C8) 
AM A EN j 


where the composition of the infalling debris is expressed by the left-hand side. In the decreasing abundance phase, 
the mass of each element in the white dwarf gradually diminishes exponentially, so that if accretion ceases at to, then 
at some time t later, the photospheric mass ratios are given by 


Xzi(t—to) ` Kalte" Hal 


= . Co 
Kat — to)  Xz2(to)e €- t0)/7z2 (C9) 


Figure 6 demonstrates that both the increasing (model 1) and steady-state abundance phases (model 2) are not 
a good match to the observational data as exemplified by Z/Mg. Instead, there is a clear and correlated pattern 
of under-abundance with atomic mass, and thus diffusion time, for both scenarios. Even for the steady state, the 
chondritic model yields x2 = 22.5, and is a poor fit. These results indicate that WD 0816-310 is observed in the 
decreasing abundance phase. 

Adopting a model where chondritic material is allowed to chemically segregate via differential sinking using Equa- 
tion C9, two cases were considered for the initial conditions X,;(to)/X,2(to): model 3 assumes an increasing phase 
at to, whereas model 4 adopts a steady state at tọ. Post-accretion, both scenarios provide excellent fits to the data, 
where x2 = 0.35 for model 3 and a decreasing phase lasting 0.76 Myr, and x2 = 0.10 for model 4 with a decreasing 
phase of 0.31 Myr. The latter model is adopted because it provides the closest agreement with the data. 


D. CONDITIONS IN ACCRETION REGION 


First, we consider conditions in the accretion region closer to the star than the transition region where solids rapidly 
sublimate, within about 5 R,, in the absence of a magnetic field. In this region, we expect a kind of cyclone, with gas 
flowing around the star in a disk with close to local orbital speed. The circular orbital velocity at distance r from the 
white dwarf is given by 

Vorb = (GM, Zell? (D10) 


With M, = 0.57Mọo and R, = 0.013Ro, the orbital velocity ranges from about 1300kms~! in the region where solids 
rapidly sublimate to about 3000 km s71 close to the star. With 1.60 x 1071? ergeV~! and metal atoms of typical atomic 
weight 35, the kinetic energies of the individual atoms are in the range of about 0.1-2.0 MeV. 
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If no magnetic field is present, one might expect the high-speed gas flow around the star to settle into a disk in 
approximate hydrostatic equilibrium normal to the plane of the disk. The ion density n;(z) as a function of distance 
z from the disk plane, where the ion density is non, would be then given approximately by 


ng(z) = nioe?’ , (D11) 


where gan is the atom plus ion number density at disk mid-plane; the scale height H is set roughly by the component 
of the stellar gravity normal to the disk, and is given by 


QkT S 


H = VOM, Amp) 


(D12) 
where k is the Boltzmann constant, T is the ambient temperature in the gas disk, G is the gravitational constant, 
A ~ 35 is the mean atomic weight of the gas atoms, Mp is the proton mass, and r is the distance from the white dwarf. 
If d is the distance from the centre of the white dwarf in units of the white dwarf radius, we find that H ~ 3d3/? km, 
which ranges from about 3km at the stellar surface to 30km in the region where solids rapidly sublimate. 


We use the equation of continuity 
M ~ nio(d) (27 Ry) (2H) Vacc(Amp) , (D13) 


where Vace is the effective radial drift velocity of atoms towards the stellar surface, to estimate the order of magnitude 
of the ion density nio(d) in the gaseous accretion radius. Assuming Vace ~ Uth, Where Gu is the mean local thermal 
velocity dispersion of atoms within the accretion disk), we find 


nio(d) ~ 5 x 10°Mgd~*/? (D14) 


1. The particle density in the 


particles per cubic centimeter, where Msg is the mass accretion rate in units of 108 gs7 
gaseous disk is thus rather low. 

If the white dwarf has a magnetic field, and the infalling gas completely un-ionized, the situation does not change. 
The gas circulates in tropical storm mode around the star. However, if an atom becomes ionized, it immediately feels 
the Lorentz force and begins to circle a local magnetic field line. The radius ry of the new orbit around the field line 
is given in Gaussian cgs units by 


Tt zs Ampvorve/qB (D15) 


where q is the charge of the electron. For a space velocity of Act, ~ 2000kms~!, and in a field of 1kG, this radius is 
of the order of 100 cm, so the ion is closely confined to a single line of force. 

A newly created ion will thus find itself moving relative to the neutral ions in the gas disk, with a relative velocity of 
the order of 108 cms™!, effectively carrying a kinetic energy relative to the neutral atoms of the order of 1 MeV. If the 
local neutral particle density is of the order of 10° cm~°, the new ion will suffer a few collisions per second. Most of 
these collisions will ionize the neutral particles struck, and these new ions will join the set locked to field lines, moving 
at high velocity relative to the neutrals. It thus appears that formation of even a few ions in the inflowing gas will 
quickly lead to a cascade of ionization, and the gas will quickly become ionized overall. In this situation, the inflow of 
gas to the stellar surface is expected to be forced to follow field lines, leading to accretion mainly around the magnetic 
poles rather than in an annulus around the orbital plane of the accretion disk. 


DATA AVAILABILITY 


FORS2 data are available through the ESO public archive (program ID: 110.243A.001, PI=Bagnulo). ISIS data are 
available at the Astronomical Data Center. TESS data used in this paper can be found in MAST: 10.17909/3es1-hp33 
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